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The two particle inclusive cross section for the reaction p + p — > 7 + c + X is 
studied in perturbative quantum chromodynamics at order O(a^). Differential 
distributions are provided for various observables, and a comparison is made with 
preliminary data from the CDF collaboration. 



1 Introduction 

The CDF Collaboration! are analyzing data on prompt photon production in 
association with charm decay products such as , jji^ or . This two-particle 
inclusive reaction is particularly interesting because it offers the possiblility 
of a detailed study of the underlying QCD dynamics such as, e.g., rapidity 
correlations. In addition the data may provide a direct measurement of the 
charm quark density of the proton due to the dominance of the eg scattering 
process. 

We completed two next-to-leading order pertuxhative QCD calculations 
of the reaction p + p^^f + c + X at high energyEru . In these calculations 
two different techniques were used in performing the phase-space integrals. In 
the first, purely analytical techniques were used. In the second approach, we 
used a combination of analytical and Monte Carlo techniques, which is more 
flexible and allows implementation of isolation cuts and other experimentally 
relevant selections. To warrant use of perturbation theory and the massless 
charm approximation, we limited our considerations to values of transverse 
momenta of the photon and charm quark p^' c > 10 GeV. In these calculations 
we consider and include hard scattering processes in which the charm quark 
is a constituent of the incident hadrons e.g., eg — > ^cX. Details of these 
calculations can be found in the above references. 
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Figure 1: Cross section as a function of z at y"' = in (a) leading order and (b) next-to- 
lcading showing the dominant contributions. 



We recently completed a different calculation of this reaction in which 
we included a finite charm quark mass explicitly. This calculation, which 
should be valid for p9p < 10 GeV, was done in leading order QCD (0(a 2 s )). 
In this case the only subprocesses contributing to the reaction are those in 
which the charm quark is produced in the hard scattering reaction, namely 
gg — > cc7 and qq — > cc-f. There are no collincar or soft singularities at this 
order and no fragmentation contributions to the reaction. The integration over 
the final state phase space could therefore be done numerically with Monte 
Carlo routines and photon isolation cuts could be implemented. 



2 Numerical Results 

Our results are presented at a center-of-mass energy yfs = 1.8 TcV. Rcnor- 
malization/factorization scales are taken as /j, = Pj,. We sum over charm and 
anticharm production throughout. Some results are presented as distributions 
in the variable z = Pt-Pt/(Pt) 2 > ^ or nn ite bin widths of z. In Fig. la the net 
lowest order contribution for the massless calculation is shown as a function of 
z at pj, = 15 GeV and for a bin size Az = 0.2. The lowest order cross section 
is made up of the lowest order direct term eg — > 7c, which is proportional to 
5(1 — z) and provides the peak at z = 1, plus the various photon fragmenta- 
tion contributions which contribute in the region z > 1. The most significant 
feature of this curve is that there is no contribution to the cross section in the 
region z < 1. This unrealistic prediction shows the inadequacy of the lowest 
order predictions. 
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Fig. lb shows the distribution in z predicted by the next-to- leading order 
massless calculation. The next-to-leading order contributions serve to lower 
the peak at z = 1, and they broaden the distribution. The cross section is 
finite at all values of z, closer to the situation observed in experiments. In 
addition, in Fig. lb we display contributions from the most important subpro- 
cesses. The eg initiated process dominates the cross section, but there are 
important contributions from the gg and cq initiated process in the IohljjJ 
region. Predicted distributions in p T , p T and rapidities may be found in Era. 

In Fig. 2a we show the z distribution for the same kinematic variables as 
above but now for the leading order massive calculation with various choices for 
the charm quark mass. For comparison the full massless next-to-leading order 
results are also displayed. Although the massive calculation is done in leading 
order, the contributing subprocesses are two-to-three scattering processes, and 
they contribute in the region z < 1. In this sense the massive calculation 
gives a more realistic result than the corresponding massless calculation. On 
the other hand it clearly predicts a smaller cross section than the next-to- 
leading order massless calculation. We note, however, that the next-to-leading 
order contributions to heavy quark production are known to provide significant 
increases in the predicted yields □. Fig. 2a also demonstrates the sensitivity to 
the choice for the charm mass particularly in the region around z = 1. This is 
the region in which the finite charm mass regulates the initial state collinear 
divergences. 

In Fig. 2b we compare our results to CDF datai for photon plus ^ pro- 
duction. The three upper points are obtained^ from the Monte Carlo event 
generator Pythia whereas the lower point is that given by our theoretical cal- 
culation. The Pythia cross sections lie substantially below the data whereas 
our cross section exceeds the data but somewhat closer to it. 



3 Conclusions 

We presented the results of three calculations of the inclusive production of a 
prompt photon in association with a heavy quark. Two analyses are done at 
next-to-leading order in perturbative QCD in the massless charm framework, 
and one calculation in leading order QCD in the massive charm framework. 
Our results in the massless cases agree quantitatively, as they should, but 
the combination of analytic and Monte Carlo methods is more versatile. A 
comparison of our next-to-leading predictions with the preliminary CDF data 
shows reasonable agreement. 
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Figure 2: (a) z distribution of the leading order massive calculation for various values of 
the charm mass. The next-to-leading order massless results are also shown, (b) Ratio of 
the measured cross section to that predicted by theory for various final state charm decay 
products. 
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